Gray- Tone Lithography Implementation of Drexhage's Method for Calibrating 
Radiative and Nonradiative Decay Constants of Fluorophores 
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We present a straightforward method to realize non-planar dielectric structures with a controlled 
height profile for use in calibration of fluorophores. Calibration of fluorescence quantum efficiency 
and intrinsic radiative and nonradiative decay rates of emitters is possible by using changes in the 
local density of optical states, provided one can control the emitter-surface distance with nanometer 
accuracy. We realize a method that is accurate yet fast to implement. We fabricate PMMA wedges 
(4mmx4mmx2 urn) by gray-tone UV-lithography of Shipley Sf8f3G2. Its applicability as dielectric 
spacer is demonstrated in Drexhage experiments for three different emitters in the visible and near- 
infrared wavelength regime. The decay-rate dependence of the fluorescent state of emitters on the 
distance to a silver mirror is observed and compared to calculations based on the local density of 
states. Quantitative values for (non) radiative decay rates and quantum efficiencies are extracted. 
Furthermore, we discuss how Drexhage experiments can help to scrutinize the validity of effective 
material parameters of metamaterials in the near field regime. 
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INTRODUCTION 

A key contribution of photonic technology to soci- 
ety is the realization of light sources with desirable 
properties, such as controlled spectrum, brightness, 
improved wall-plug efficiency, or coherence proper- 
ties. The creation of such novel light sources, i.e., 
LED's and lasers, heavily depends on two types of 
innovation. These are firstly the realization of cheap 
materials that have optical transitions of high effi- 
ciency and desired frequency, and secondly the in- 
tegration of these source materials with photonic 
structures that further improve the performance. 
Among the material developments, innovations in 
the last decade range from organic dyes and light- 
emitting polymers in organic LEDs to the use 
of epitaxially grown III-V semiconductor quantum 
dots 0] and quantum wells 0] in common laser 
diodes and LED's, to colloidally fabricated II- VI 
quantum dots [H that are size tunable throughout 
the entire visible spectrum. To improve the per- 
formance of the bare sources, or homogeneous thin 
layers of the luminescent materials, many groups 
pursue nanophotonic techniques to improve light 
extraction, as well as to increase the spontaneous 
emission rate in favor of undesired nonradiative de- 
cay constants. Both, photonic crystals [sj as well as 
patterning with arrays of plasmonic antennas @-@] 
have been shown to succesfully contribute to light 
extraction and optimization of the radiative time 
constants [10|. The key quantity that quantifies the 
possible improvement of the radiative decay is the 
local density of optical states (LDOS) enhancement 
that a nanophotonic structure provides [O, E3] • The 
LDOS counts the total number of optical modes, 



weighted by their amplitude right at the emitter, 
that are available to the source for radiative transi- 
tions, and hence directly appears in Fermi's Golden 
Rule for spontaneous emission . Optimization of 
LDOS is not only important for macroscopic classi- 
cal light sources, such as LEDs, but is also key for 
quantum optical applications where a single quan- 
tum emitter is to be strongly coupled to a single 
optical mode (13 . [13] • 

A method to reliably measure the intrinsic time 
constants of arbitrary emitters is of large importance 
both for quantifying improvements in emissive ma- 
terials, and for quantifying the LDOS enhancing po- 
tential of a nanophotonic structure. For an emissive 
material, the challenge is to determine rapidly and 
accurately the intrinsic radiative rate, nonradiative 
rate and quantum efficiency of emitters as a tool to 
guide material improvement and to understand the 
mechanisms behind, e.g., unwanted nonradiative de- 
cay. Conversely, if one wishes to benchmark the lo- 
cal density of states improvement that a structure 
can intrinsically provide, it is important to probe the 
structure with a source that has first been quanti- 
fied in terms of its radiative and nonradiative rate 
constants. Unfortunately, commonly used fluores- 
cence decay measurements only provide the sum of 
rates [l~5j |. while quantum efficiency measurements 
are usually based on brightness comparisons. Such 
comparisons are prone to imprecision if one has to 
rely on comparison to a fluorescence quantum effi- 
ciency standard, and yield erroneous results when 
a sample consists of a heterogeneous ensemble in 
which fluorophores exhibit large brightness varia- 
tions. On basis of an experiment first performed by 
Drexhage [lf^ [13] > many authors have realized that 
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intrinsic rate constants can be reliably measured by 
applyin g a known, controlled LDOS variation to an 
emitter [li-il. Drexhage studied the radiative de- 
cay rate of europium ions as a function of distance 
to a silver mirror. The observed variation in decay 
rate, stemming from interference of emitted and re- 
flected light, can be explained by a change in the 
LDOS at the emitter position [2l|. Since the LDOS 
is exactly known, the data can be quantitatively 
separated into an intrinsic nonradiative rate that 
does not vary with LDOS, and a radiative rate that 
does. This technique has been used for organic dyes, 
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rare earth ions 113 . 1221 ] , and more recently for II- VI 
quantum dots j23l ]. III-V quantum dots 2J, 25], and 
even single emitters when using a nanomechanically 
scanned mirror [23]. Unfortunately, the techniques 
to controllably vary distance are generally elaborate, 
and material specific. For instance, Leistikow et 
al. (23[ | required fabrication of a large set of sam- 
ples with evaporated layers of calibrated heights. In 
the case of Stobbe et al. [24|, a single substrate could 
be used, but an elaborate reactive ion etching step 
specific for III-V chemistry, and using a complicated 
masking procedure was required to fabricate discrete 
steps. In this work we report an easily implemented 
method to realize Drexhage experiments on top of 
arbitrary planar structures. 

In this paper we propose that gray-tone opti- 
cal lithography [23] allows to attach very shallow 
wedges on top of arbitrary substrates. Drexhage 
experiments can then be performed either by de- 
positing the wedge on the emitter and evaporat- 
ing a mirror on it, or conversely by placing the 
wedge on a mirror, and distributing sources on the 
wedge (Fig. [T|). The key requirement for the opti- 
cal wavelength regime is that the wedge has a shal- 
low angle so that the mirror is almost parallel to 
the emissive substrate, yet also to have nanome- 
ter control over the wedge height and roughness 
that sets the spatial separation of the emitter and 
the substrate. We fabricated wedge-shaped dielec- 
tric spacers by gray-tone lithography and performed 
Drexhage experiments to calibrate three different 
emitters: fluorescent polystyrene beads emitting 
at 605 nm, CdSeTe/ZnS (core/shell) quantum dots 
emitting at 800 nm, and Dibenzoterrylene molecules 
in anthracene crystals emitting at 750 nm. While 
these emitters were chosen for their promising use 
in nanoscale quantum optics with plasmon anten- 
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- 301 ] and metamaterials [3 114331] . the method 
is easily applied to any emitter that can be reason- 
ably homogeneously distributed in a planar layer. 
The wedges can also be used with calibrated emit- 
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FIG. 1. A wedge shaped dielectric separates a layer of emit- 
ters from a mirror. This geometry allows to conduct a Drex- 
hage experiment — observing the fluorescent decay rate of 
emitters as a function of distance to a mirror — on a single 
sample. 



ters, to measure LDOS near unknown substrates, 
such as gratings and plasmon antenna arrays. The 
paper is structured as follows. In the fabrication 
section we discuss gray-tone UV-lithography, wedge- 
profile retrievement and emitter layer spin-coating 
steps. The experiment section covers the optical 
measurement procedure. Lifetime measurements 
and derived intrinsic rate constants for the three 
emitters are presented in the results section. Finally, 
along with our conclusion we provide an outlook. 

FABRICATION 

In this paper, we discuss four samples: three sam- 
ples featuring emitters with specific emission wave- 
lengths and fluorescent lifetime characteristics and 
one benchmark sample to check the height pro- 
file using a variety of techniques. All samples fea- 
ture an optically thick (~ 100 nm) layer of silver 
that is evaporated on silicon wafer pieces (each 
~ 20 mm x 20 mm). The silver layer, which we char- 
acterized by ellipsometry, serves as mirror for the 
Drexhage experiments. The essential step is the 
fabrication of a dielectric spacer with controlled 
and graded height profile using optical lithography. 
Since optical lithography resists tend to fluoresce 
themselves, we create two types of samples, namely 
samples with and without an intermediate, consid- 
erably less fluorescent, PMMA layer. For samples 
of the first type a 2 um thick layer of PMMA (M = 
950 000g/mol, 8% in Anisole) is spincoated (45 s at 
1000 rpm, baking for 2min at 180°C). For all sam- 
ples (with and without PMMA) we then spincoat a 
2 |j.m thick layer of Shipley Microposit S1813 G2 pos- 
itive UV-resist (45 s at 1000 rpm, baking for 2 min at 
115°C). To define the wedge shape in the resist, we 
perform UV-lithography using a Stiss MJB3 mask 
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FIG. 2. Diffusive UV-illumination of a 'gray-tone' chromium mask generates a wedge shaped zone with doses above the critical 
dose of the positive resist. After development, the wedge can be transferred into a different material, such as PMMA, via 
reactive oxygen ion etching (a). Differential interference contrast microscopy images of developed Shipley S1813G2 resist 
without diffusive element during UV-exposure, resembling the chromium mask pattern (b), and with diffusive element, showing 
a 'smoothed out' continuously sloping surface (c). The interference color change along x implies varying distances to the silver 
mirror. 



aligner with a binary chromium mask that consists 
of parallel lines, at varying surface coverage similar 



to Christophersen et al. [26f] . The mask for each 



dielectric wedge is a 4 mm x 4 mm area, made up of 
4 mm long chromium lines in the y-direction with in- 
creasing width in the x-direction from 1.5 (J.m (low 
end of the wedge) up to 12 um (high end of the 
wedge). All lines are spaced in the x-direction by 
gaps of 3 urn. Thereby, the average density con- 
tinuously varies from 33 % surface coverage to 80 % 
surface coverage. To generate a graded illumina- 
tion, we use an opal glass diffuser (Edmund Optics, 
NT02-149) placed in a filter mold ~ 1 cm above the 
sample and mask. The resist is exposed with a dose 
of 300mJ/cm 2 and developed in Microposit MF-319 
for 15 s (Fig. [2]). At this point, we have a wedge pro- 
file in S1813 on top of Ag or a Ag/PMMA stack. 
For the two samples with PMMA, the wedge pro- 
file, defined in S1813, is transferred into the PMMA 
layer by reactive oxygen ion etching (Oxford Instru- 
ments Plasmalab 80+, 20 seem 2 gas flow, 50 W 
forward power, 292 K operating temperature). Af- 



ter approx. 35 min of etching the S1813 is completely 
removed from the lower end of the wedge. The etch 
rate of PMMA is twice the etch rate of S1813, there- 
fore the slope of the original profile is changed and 
the upper end of the PMMA wedge will be covered 
with residual S1813. 

We characterize the wedge profiles by profilome- 
try (KLA-Tencor Alpha-Step 500) and compare the 
profilometry data set of the benchmark sample to 
Mirau-interferometry data (Nikon 20x CF IC Epi 
Plan DI) (Fig.[3|a)). Since the lateral extent of the 
wedges is large (4 mm) , raw profilometry data shows 
curvature and roughness that are due to the scanner, 
rather than the sample. However, once one uses a 
flat optical substrate, such as an optical grade silver 
mirror (Thorlabs PF10-03-P01), as reference in the 
profilometer, a good agreement between profilome- 
try and interferometry is obtained (Fig.[3|a), dots 
and curve). This supports the sole use of the pro- 
filometer as a tool to acquire (x,z)-profiles of all the 
other samples. Typically, the base of a wedge, in 
contact to the silver mirror, is 4 mm wide in the 
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FIG. 3. Panel (a): height profiles of 4 different wedges, used for fluorescent beads, quantum dots, DBT, and a topographic 
benchmark sample. For the benchmark we also performed Mirau-interferometry (round symbols) and find good agreement 
with profilometry. The fluorescent bead and DBT in AC samples feature an additional PMMA layer underneath the S1813 
UV-resist. The wedge for the quantum dot sample consists solely of S1813 (a). Panel (b): an atomic force microscopy scan 
is conducted to retrieve a roughness estimate of the typical wedge surface. Line traces in (c) show a sub-30 nm roughness as 
indicated by the double-arrow. 



x- and y- (scan-) directions, as inherited from the 
mask. The (x,z)-profiles resemble a 'fin'-shape start- 
ing off with a steeper slope at the low end and end- 
ing with a shallower height increment per sideways 
displacement at the upper end of the wedges. The 
overall slope remains well below 1.5 (a.m/mm, mean- 
ing that the wedge angle is so shallow that effectively 
constant-height data can be obtained in micro- 
fluorescence experiments. Since neither profilom- 
etry nor Mirau-interferometry can resolve rough- 
ness reliably on lateral length scales below 1 (a.m, 
we also performed atomic force microscopy mea- 
surements. We extract a typical surface roughness 
of Az/Ax < 30nm/50um (Fig.[3(b),(c)). Com- 
pared to established methods for Drexhage exper- 
iments, which are usually based on controlled step- 
wise RIE etching, or on controlled evaporation of 
different thicknesses, our method has the advantage 
that all desired source-mirror distances are created 
on a single substrate with a single simple optical 
lithography step. The roughness and shallow slope 
of the wedge imply that Drexhage experiments can 
be performed for wavelengths down to ~ 500 nm 
(< A/10 roughness). Finally, we note that Fig. [3(a) 
shows that while all samples have a smooth wedge, 
not all wedges have the same overall height profile. 



In general etching from S1813 into PMMA steepens 
the height profile due to a factor two difference in 
etch rate. In our work differences in height profile 
also occur due to minor variations in the spincoated 
thickness and in the UV exposure dose, as a re- 
sult of adjustments in the mask aligner illumination 
alignment between runs. In this work we use indi- 
vidual z-profile calibrations for each wedge, though 
improvements in processing could render this proce- 
dure unnecessary. 

A succesful Drexhage experiment not only re- 
quires smooth dielectric spacers of controlled height, 
as characterized above, but also that the emit- 
ters can be homogeneously dispersed on the wedge. 
We present measurements on three types of fluo- 
rophores. Firstly, we use fluorescent polystyrene 
beads (Invitrogen Fluospheres F8801), with a nomi- 
nal bead diameter of 0.1 um containing about 10 3 
randomly oriented dye-molecules each, that are 
completely chemically shielded from the bead's envi- 
ronment. We estimate this number of dye molecules 
on the basis of intensity measurements taken in a 
single molecule sensitive microscope of known col- 
lection efficiency as reported by Primmer et al. [Hjl . 
Their fluorescence intensity peaks at 605 nm. Since 
S1813 wedges in themselves fluoresce in the same 
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wavelength range, with comparable time constants, 
we use a nonfluorescent PMMA wedge. To coun- 
teract agglomeration of beads we sonicate the stock 
solution (2 % solids) for 2 minutes prior to mixing 

1 ial of the bead solution with 1 ml of deionized wa- 
ter. This mixture is spin-coated for 10 s at 500 rpm 
(lOOrpm/s acceleration) followed by a second spin 
step of 45 s at 1500 rpm (500 rpm/s). Prom a typical 
confocal fluorescence microscopy image (Fig.[5l left) 
we deduce a bead concentration of ~ 0.15 um~ 2 . 

As a second emitter, we study CdSeTe/ZnS 
(core/shell) quantum dots (Invitrogen Qdot 800 
ITK carboxyl) as they were recently proposed as 
light sources for plasmonic applications [23, [2§] . In 
comparison to more common CdSe/ZnS quantum 
dot nanocrystals for visible emission, these quantum 
dots feature longer lifetimes [35[ . The quantum dots 
are diluted in a borate buffer solution (pH = 8.0) to 
a molar concentration of 8 (xM. We mix 4 \il of this 
solution with 1 ml of deionized water and spincoat 
a droplet of the mixture on top of the S1813 wedge- 
sample for 45s at 4000 rpm (4000 rpm/s). Here, the 
use of an S1813 wedge is justified, since the fluo- 
rescence of the quantum dots and of the polymer 
can be separated easily via spectral selection, and a 
factor 10 contrast in fluorescence decay constants. 

As third type of emitter, we investigate Diben- 
zoterrylene (DBT) molecules in anthracene (AC) 
crystals [3^, [37j • DBT molecules have been found 
to feature great photostability and brightness at 
room temperature (10 12 emitted near-infrared pho- 
tons before photobleaching) while also fluorescing 
in the near-infrared around 750 nm [37j. We dis- 
solve 0.6 mg of DBT powder (Dr. W. Schmidt, 
PAH Reasearch Institute, Greifenberg, Germany) 
in 1 ml of Toluene. In accordance with the recipe 
by Toninelli [37j |. we use 5.3 mg of AC dissolved in 

2 ml of Diethylether and 20 \x\ of Benzene. Finally, 
80 ul of DBT-solution is mixed with 0.5 ml of AC- 
solution. The resulting solution is spin-coated for 
40s at 3500 rpm (500 rpm/s acceleration). 



EXPERIMENT 

Fig. d] depicts the basic components of the con- 
focal fluorescence lifetime scanning microscope, in 
which we perform the Drexhage experiments. Light 
from the excitation laser is focused to the diffrac- 
tion limit, and emitted light from fluorescent ob- 
jects in the focus is collected via the same objec- 
tive (Nikon lOOx, NA = 0.90, Plan Fluor) in a con- 
focal arrangement. The fluorescent light is sepa- 
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FIG. 4. Sketch of the time correlated single photon counting 
setup. Light from a pulsed pump laser is focused onto the 
emitter layer. A dichroic mirror is used to separate the fluo- 
rescent light stemming from the emitters and the pump light. 
Depending on the emission wavelength of the emitters, appro- 
priate long-/bandpass filters are added to the beam path. In 
that way only light stemming from fluorescent objects of in- 
terest is focused onto the APD. The photon arrival times with 
respect to the pump pulse together with the (x,y)-piezo stage 
position generate a fluorescent lifetime image of the sample 
at a certain emitter-mirror separation z. 



rated from the excitation at a dichroic beamsplit- 
ter and passes additional long- /bandpass filters cho- 
sen according to the absorption and emission spec- 
trum of the emitter of interest. The fluorescent 
beads as well as the quantum dots are pumped by a 
pulsed laser (Time-Bandwidth Products) operating 
at 532 nm emission wavelength (green), 10 MHz rep- 
etition rate with pulse widths < 10 ps. Dibenzoter- 
rylene molecules have a 25 times higher absorption 
cross section in the red compared to the green part 
of the spectrum. Therefore, we choose a different 
pump source for this emitter: a pulsed laser diode 
(Edinburgh Instruments EPL) at 635 nm, operated 
at 10 MHz repetition rate featuring pulse widths of 
< 100 ps. 

Achromatic optics focuses the fluorescent light 
onto a silicon avalanche photodiode (APD) (ID 
Quantique idl00-20ULN) with an active area di- 
ameter of 20 (j.m, using an effective magnification 
from objective to APD of 20x. The APD-pulses 
and the reference pulses from a trigger diode (green 
laser) or electrical trigger output (red laser) are reg- 
istered by a picosecond, 16 channel, pulse correlator 
(Becker&Hickl DPC-230), which records absolute 
timestamps at 165 ps resolution for each event. The 
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FIG. 5. Typical fluorescence intensity image of beads on top 
of the PMMA wedge (left) and corresponding fluorescence 
decay histogram (accumulated events from diffraction limited 
bright spots in intensity image) (right). We clearly see a sin- 
gle exponential decay for which intensity and lifetime r are 
fitted in the maximum likelihood sense assuming Poissonian 
counting statistics for each time bin. 



sample is mounted on a (x,y)-piezo stage. By scan- 
ning the sample with respect to the objective we ac- 
quire fluorescence intensity maps of 100 x 100 pixels 
(~ 10 (a.mxl0 (o.m) with a pixel scan rate of 100 Hz. 
For each scan we sum the single photon events of a 
chosen region of interest in a time-histogram. These 
areas are so small that no appreciable height gradi- 
ent occurs in the wedge. The images hence serve 
to assess local homogeneity in fluorophore intensity 
only. For each scan we correlate the single photon 
detection events and laser pulses to form a fluores- 
cence decay histogram, which we sum either over 
the full image or over a chosen region of interest 
(e.g., to select isolated beads) in the 2D scan. Differ- 
ent positions along the height gradient of the wedge 
are reached by a manually operated micromechani- 
cal stage. 



RESULTS 

Fig. left shows a typical scan of beads on top of 
the wedge aquired by fluorescence-lifetime imaging 
microscopy (FLIM). The intensity map shows iso- 
lated diffraction limited bright spots on a very dark 
background, which is consistent with having a dilute 
sprinkling of isolated beads. Although the slight 
top-bottom intensity gradient in this particular im- 
age suggests a small focus drift during the mea- 
surement, our method is robust against such drift. 
Raster scanning excludes photophysical changes in 
the sample, such as bleaching as artifact. Since our 
quantum efficiency calibration does not rely on ex- 
tracting brightness, but on extracting lifetimes, we 
are not prone to, e.g., the small setup drifts that 
can cause brightness artifacts. Selecting only those 
events that appear in the diffraction limited bright 



spots we obtain a histogram of photon counts versus 
arrival time after the pump pulse, i.e., a fluorescence 
decay curve, for all beads within the scanned image. 
One such histogram is shown in Fig.[5l right, taken 
from the fluorescence image Fig.0 left recorded 
around position x = 340 urn along the wedge length 
(i.e., z ~ 700 nm). For each position along the 
wedge profile, i.e., for each sample-source distance 
where we acquired a FLIM image, we fit the ob- 
served fluoresence decay histogram with a single ex- 
ponential decay. In our fit we assumed the back- 
ground to be fixed (set to the mean events per bin 
in the time bins with time delay < 0, before pump 
pulse arrival) and fit only a decay rate and an initial 
amplitude. We find the most probable set of these 
two parameters by using the maximum-likelihood 
method under the assumption of Poissonian count- 
ing statistics for each time bin [38]. In Fig.[6l left we 
plot the experimentally retrieved fluorescent decay 
rate as a function of bead-mirror distance z. We 
find an increase of the decay rate up to 0.27 ns -1 
for z « 90 nm. We did not resolve decay rates 
for values closer to the mirror which, predicted by 
LDOS theory, would first decrease to a minimum at 
around z = 25 nm before rising again for z — > 0. 
For z > 100 nm, the decay rate dependancy resem- 
bles a damped oscillation around 0.18 ns -1 with a 
periodicity of « A£™ ad /(2n PM MA) ~ 200 nm as ex- 
pected when considering the interference condition 
for emitted and reflected field amplitudes. 

In order to extract a quantitative calibration of 
the intrinsic decay constants of the emitters, i.e., 
their intrinsic nonradiative decay rate and their ra- 
diative decay rate when held in vacuum, we fit the 
experimentally acquired total decay rate ^ to t(z), de- 
pendent on the emitter-mirror separation z, to 



ltot(z) = 7rad ■ P{z) + 7nonrad 



(1) 



where p{z) equals the relative local density of states 
p =LDOS/LDOS vac for a dipole above a silver mir- 
ror. Note that the intrinsic nonradiative decay 
7nonrad acts as an offset that is not affected by the 
mirror. Any quenching induced by the mirror is con- 
tained in the term proportional to the intrinsic ra- 
diative decay rate, i.e., in p. We use established 
theory 12 , 1|| 2l|, [3{| to calculate p on basis of the 
refractive indices for S1813, PMMA and the silver 
mirror that were acquired from ellipsometry mea- 
surements. 

Since in our experiment, we acquire fluorescence 
intensity from an ensemble of ~ 10 3 dye molecules 
within each bead, we assume an isotropic distribu- 
tion of dipole orientations. Therefore we use the 
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orientation averaged 



Piso(z) = \p±{z) + rPuW 



(2) 



where p±(z) and p±(z) are the relative local den- 
sity of states for dipoles oriented perpendicular and 
parallel to the mirror surface, respectively pZl ], 

The fitted Y t s ° t (z) (calculated for dipoles embed- 
ded 10 nm below the PMMA/air interface to ac- 
count for the bead material) is plotted as the red 
line in Fig.[6l left. We find fair agreement with the 
experimental data. For distances exceeding the vac- 
uum emission wavelength of the emitter (due to the 
wedge index this encompasses several oscillations in 
the LDOS) the oscillation amplitude of the local 
density of states decreases and becomes compara- 
ble to the uncertainty in fitted total decay rates. 
We hence use the emitter's emission wavelength as 
upper bound on distances plotted. By plotting de- 
cay rate versus LDOS, instead of decay rate versus 
z, we can directly find the radiative and nonradia- 
tive decay rate from the slope and ordinate intersec- 
tion of the linear dependence, respectively (Fig. [7]). 
Throughout this work, error bars on rates and ef- 
ficiencies result from the LDOS fit, taking into ac- 
count error bars on the decay rate fitted at each 
vertical distance. For the fluorescent beads, we find 
an intrinsic radiative decay rate of (0.10 ±0.01) ns _1 
and a nonradiative decay rate of (0.07 ± 0.01) ns _1 . 
The hypothetical radiative decay rate for this emit- 
ter in vacuum, extrapolates to a radiative rate of 
(0.15 ±0.01)ns~ 1 in its PMMA host. We extract a 
quantum efficiency 



Q'&' — Trad/ vTrad ± 'Ynonrad) 



(3) 



of 61% ± 7% for the emitter hypothetically in vac- 
uum, and around 70% ± 6% for the emitter embed- 
ded in its bulk host material, i.e., bulk PMMA. 

In many applications targeting the use of meta- 
materials and plasmonics to control emitters, one 
preferentially does not use emitters with emission 
wavelengths at 605 nm, such as the beads, but rather 
emitters emitting further into the near-infrared. Re- 
cent reports by Curto et al. (2(| propose that Cd- 
SeTe/ZnS quantum dots emitting at around 800 nm 
are ideally suited single emitters for plasmonic ap- 
plications, as they are ultrabright. The fluorescence 
time delay histogram retrieved from a FLIM image 
of an ensemble of Invitrogen Qdot 800 ITK carboxyl 
quantum dots does not resemble a single exponen- 
tial decay. Here, we make use of the fact that an 
ensemble of quantum dots should be modeled with a 
continuous distribution of decay rates 4l|, |42J . More 



precisely, the natural logarithm of decay rates 7 is 
assumed to be normally distributed according to the 
log-normal distribution 



p(-y) = A ■ exp 



In 2 (7/7" 



(4) 



The normalization constant A is given by the condi- 
tion J °°p(7)d7 = 1. The dimensionless width w can 
be rewritten as the width of the rate distribution for 
which p = 1/e: 



A7 = 27 mf sinhu;. 



(5) 



Hence, the two free parameters of our fit-model are 
the most frequent decay rate 7 mf , at which the log- 
normal distribution is centered, and the rate distri- 
bution width A7. The fitted most frequent decay 
rates are shown in Fig.(6j center and Fig. cen- 
ter together with a fit to the LDOS for isotropi- 
cally oriented dipoles. We would like to point out 
that as we fit the Drexhage model to just the re- 
trieved 7 mf , the extracted values we quote quan- 
tify only the most frequently occurring decay rates 
and quantum efficiencies, and not the width of the 
underlying distribution. For the quantum dots we 
find a most frequent intrinsic radiative decay rate 
of (5.46 ± 0.51) us -1 , a most frequent nonradiative 
decay rate of (0.62 ±0.66) i^s -1 and a most frequent 
quantum efficiency of 90% ±11% for the emitter hy- 
pothetically in vacuum. In its actual bulk host ma- 
terial, i.e., bulk S1813 we find a most frequent radia- 
tive decay rate of (9.00 ± 0.84) ps -1 and a quantum 
efficiency of around 94% ± 7%. 

The high quantum efficiency for these quantum 
dots makes them exceptionally promising for plas- 
mon quantum optics, since operating at 800 nm op- 
timizes the emission frequency to be at the inter- 
section of low plasmon loss, yet efficient silicon de- 
tection. The high quantum efficiency is surprising 
given that the solvent /ligand exchange to aqueous 
condition, and operation of quantum dots unpro- 
tected against oxygen as in our experiment usu- 
ally adversely affects the photophysical properties. 
Here we note that measuring quantum efficiency 
via modulation of radiative lifetime selectively mea- 
sures the quantum efficiency of the ensemble of 
dots that radiate, i.e., that are not temporarily 
or permanently dark due to (photo) chemical pro- 
cesses such as bleaching or blinking. This prop- 
erty should be contrasted to absorption/emission 
brightness measurements that report quantum ef- 
ficiency as the fraction of absorbed photons that 
is converted into radiated photons by an ensem- 
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ble of nanocrystals. In such measurements a mix- 
ture of a dark subensemble that absorbs but does 
not emit with a subensemble of unit-efficiency emit- 
ters is indistinguishable from a homogeneous sam- 
ple of emitters with below-unit quantum efficiency. 
The conclusion from our measurement is that those 
dots that are not (photo) chemically altered by an 
oxygen-rich or aqueous environment, retain high 
quantum efficiency. A similar conclusion that the 
quantum efficiency of those dots that radiate is much 
higher than the ensemble efficiency obtained from a 
absorption/emission brightness measurements was 
also reached by Leistikow [2j|, for CdSe quantum 
dots emitting around 600 nm, with efficiencies be- 
tween 66% and 89%. Evidently, the CdSeTe/ZnS 
core shell dots emitting at 800 nm have a 6- fold 
longer lifetime compared to CdSe quantum dots, 
consistent with observations by Vion et al. 35J]. It 
is remarkable that the (bright) 800 nm dots in our 
work manage to retain a very high quantum effi- 
ciency, despite this 6- fold longer lifetime [231 ]. 

An interesting alternative to the quantum dots, 
which while efficient, have a slow and non-single ex- 
ponential decay, could be DBT with an emission 
wavelength of ~ 750 nm j36l . 37], DBT was re- 
cently reported to be ultrastable and ultrabright 
as an emitter at room temperature [37j]. In con- 
trast to Toninelli's observations on single molecules, 
we observe a non-single exponential decay trace for 
ensembles of DBT molecules. Therefore, we apply 
the same analysis scheme as for the quantum dots: 
log- normally distributed decay rates [4l[. In this 
manner we find a most frequent intrinsic radiative 
decay rate of (0.05 ± 0.01) ns _1 , a most frequent 
nonradiative decay rate of (0.27 ± 0.01) ns _1 and 
a most frequent quantum efficiency of 16% ± 3%, 
when calculating with the radiative rate extrapo- 
lated to vacuum. Decay rate data and LDOS fit 
for DBT are shown in Fig.[6l right and Fig.d right. 
Since DBT is always used in anthracene, it is use- 
ful to extract the quantum efficiency in bulk an- 
thracene. Correcting for the index of refraction of 
anthracene, we find a most frequent radiative decay 
rate of (0.09 ± 0.02) ns -1 and a quantum efficiency 
around 24% ± 4%. 

The moderate to low quantum efficiency that our 
measurement retrieves for DBT is surprising, since 
Toninelli et al. [37j] have reported very high count 
rates from single DBT molecules in anthracene at 
room temperature. In saturation, they reported de- 
tecting photons at a rate of up to 0.5% of the in- 
verse lifetime. A quantum efficiency of around 25% 
is within, but at the low end of, the range of values 



that are consistent with the observation of Toninelli 
et al. (37| . depending on the actual microscope col- 
lection efficiency (state of the art: a few percent at 
800 nm) . We note that the observations of Toninelli 



et al. [37J] are strictly for individual molecules se- 
lected to be ideal in the sense of being long lived, 
which likely selects molecules from those 10% best 
incorporated in the anthracene matrix. While agree- 
ment with LDOS theory for isotropic dipole orien- 
tations is good, we found poor agreement of decay 
rate data to LDOS lineshapes for specific dipole ori- 
entations, despite reports of preferential orientation 
for this system [33]. This further highlights that 
we probe a heterogeneously distributed ensemble of 
emitters, as opposed to selecting particular emitters 
as in single molecule experiments [33]. From ob- 
servation of the film quality, we note that it was 
difficult to obtain homogeneous anthracene crystal 
films throughout the entire wedge length, both on 
the wedge material as well as on clean cover slips. 

Tab. [I summarizes the retrieved radiative and 
nonradiative decay constants and respective quan- 
tum efficiencies for all three fluorophores. 



CONCLUSIONS AND OUTLOOK 

In conclusion, we have demonstrated that gray- 
tone UV-lithography [26| provides a facile method 
to create samples to calibrate ensembles of emit- 
ters. The essential steps of this method are (1) 
gray-tone lithography to create S1813 or PMMA 
polymer wedges on flat reflective substrates, such 
as an Ag or Au mirror, (2) a homogeneous dispersal 
method to distribute fluorophores on the wedge, and 
(3) lifetime measurements along the length of the 
wedge to effectuate Drexhage experiments []j| . Pre- 
vious methods rather used controlled stepped etch- 
ing into substrates that had the light sources em- 
bedded in them, i.e., a material-specific technique 
suited for III-V sources [IH , or depended on the cre- 
ation of a large but discrete set of substrates with 
different deposition thicknesses of dielectric spacer 
layer [§3|]. Instead, our method allows to create 
a continuous wedge on a single substrate. Our 
method can thus be of large use for calibrating a 
range of fluorophores for applications ranging from 
quantum optics, to organic and inorganic light emit- 
ting diodes. We note that a large range of varia- 
tions is easily implemented. For instance, if surface 
chemistry means it is advantageous to first deposit 
the fluorophores on glass, then create a wedge, and 
then deposit a mirror, this is equally easily imple- 
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FIG. 6. Fitted total decay rates 7^t vs distance to the mirror z (black dots). White circles are fitted values to decay traces 
gathered right on top of the mirror (next to the wedge, z ~ 0) and plotted at a fixed height offset from zero for good visibility. 
The red lines indicate fitted relative LDOS pi so for isotropically oriented dipoles. 




FIG. 7. Total decay rates 7™t vs relative LDOS pi so for isotropically oriented dipoles. The red lines indicate fits according to 
eqffl 



mented. Beyond calibration of unknown emitters 
using a known LDOS, the same method is also of 
large interest to do the reverse, i.e., to measure an 
unknown LDOS using a calibrated emitter. Indeed, 
in plasmonics, metamaterials, and the new field of 
metasurfaces, one frequently encounters questions 
that revolve around the LDOS at controlled dis- 
tance. Consider for instance the optimization of 
light emitting diodes by plasmon particle array sur- 
faces (jl, HI, 0] or extraordinary transmission grat- 
ings [431 ] . The question at which distance one should 
optimally place the emitters to both enhance out- 
coupling, radiative rate, and quantum efficiency is 
of key importance, yet difficult to address experi- 
mentally. Likewise, the fundamental study of effec- 
tive medium parameters of metamaterials come to 
mind. Photonic metamaterials are artificial mate- 
rials with periodic arrangements of subwavelength 
scatterers aimed at arbitrary control of permittivity 
e and permeability /i to realize perfect lenses and 
invisibility cloaks via transformation optics |44H471] . 
For a wide range of metamaterials e and [i have been 
reported [48145 0] . showing that it is indeed possi- 
ble to obtain, e.g., effectively negative \i. Effective 
medium constants are commonly retrieved from far 



field experiments. Since metamaterials have intrin- 
sic near field benefits, it is interesting to test whether 
these effective material parameters retrieved from 
far field measures such as transmittance and re- 
flectance, are still valid in the near field. Using our 
wedge technique one can continuously sweep source 
height through the near field zone to examine the 
transition from resolving individual building blocks 



to resolving just the effective parameters [51|, |5 



The polymer wedges that we fabricate can indeed 
also be readily made on patterned surfaces, thereby 
opening the road to controllably vary the near-field 
spacing and measure the metamaterial LDOS. This 
method is complementary to elaborate near-field 
scanning methods [34]]. The loss of lateral resolu- 
tion in our method compared to near-field scanning 
is offset by the ease of use of the gray-tone lithogra- 
phy wedge technique, and the fact that in many ap- 
plications (e.g., plasmon enhanced LEDs) the only 
useful quantity would anyway be of ensemble aver- 
aged nature. 

We thank Martin Frimmer for initial help with 
lifetime measurements, Gijs Vollenbroek for assis- 
tance with sample fabrication, and Dr. Toninelli 
for kindly sharing her experience with depositing 
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Invitrogcn Fluospheres F8801 


Invitrogen Qdot 800 ITK carboxyl 


DBT in AC 


<v mf 

Trad. host 


(0.15 ± 0.01) ns^ 1 


(9.00 ±0.84) ns^ 1 


(0.09 ±0.02) ns" 1 


mf 

Tnonrad 


(0.07 ± 0.01) ns" 1 


(0.62 ± 0.66) us" 1 


(0.27 ± 0.01) ns" 1 


n p mf 
y- e -host 


71% ± 6% 


94% ± 7% 


24% ± 4% 



TABLE I. Table of emitter properties retrieved from our Drexhage experiments. As we assumed a log-normal distribution 
of decay rates for the probed Invitrogen Qdot and DBT ensembles, the stated decay rates and quantum efficiencies are the 
ones which are most-frequent (superscript 'mf'). Radiative decay rates and quantum efficiencies are quoted for the emitters 
embedded in their respective host medium: PMMA (Fluospheres), S1813 (quantum dots), and anthracene (DBT). 
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